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ABSTRACT 

We present results of Gemini spectroscopy and Hubble Space Telescope imaging of the 
3C 381 radio galaxy. Possible ionising mechanisms for the Extended Emission-Line Re- 
gion were studied through state-of-the-art diagnostic analysis employing line-ratios. 
Photoionisation from the central engine as well as mixed-medium photoionisation mod- 
els fail in reproducing both the strengths and the behaviour of the highest-excitation 
lines, such as [Ne v]A3424, He n, and [O iii]A5007, which are measured at very large 
distances from the AGN. Shock-ionisation models provide a better fit to the observa- 
tion. Expanding shocks with velocities higher than 500 km s _1 are capable of reaching 
the observed intensity ratios for lines with different ionisation states and excitation 
degrees. This model also provide a direct explanation of the mechanical energy input 
needed to explain the high-velocity line-splitting observed in the velocity field. 
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1 INTRODUCTION 

One of the most striking features found in radio glaxies (RG) has been the presence of large regions of i onised gas not 
only confined to central regions but also extending from tens to hundreds of parsecs from the host galaxies ([Fosburvl 1 19891 : 

These Ex- 



McCarthv et al. 19931 : Binette. Wilson fc Storchi-Bergamann 19961 : Solorzano-Inarrea. Tadhunter fc Axon 2001 



tended Emission-Line Regions (EELR £| received their name for be ing located tow ards the outskirts of the galaxies, and far 



beyond the Narrow-Line Region (NLR), which i s closer to the A GN IjPeterson 1997). However, while the NLR has shown itself 
as a direct consequence of the nuclear activity (|Petersonlll997l . and references therein), the or igin of the EELR's gas and the 
mechanisms that trigger its emission are still subject of debate ( Solorzano-Inarrea et al. 2002f ). 

These EELRs are characterised by emission of narrow forbidden lines of several species with a variety of ionisa tion states 
and excitation degrees, together with strong Hydrogen recombination lines I McCarthy et al. 19931 : Peterson 1997 ). Also, the 
shapes, intensities and morphology of this gas vary from object to object. 

There al ready exist stud i es tha t could give some clues to the physics involved in the emission of these large regions. 
For example. iMcCarthv et al.l (1 19871 ) : IChambers. Milev fc van Breugell (11987) showed that the UV continua of powerful radio 
sources at z > 0.6 are closely aligned with their radio sources. For lower redshifts and l ower luminosity radio ga laxies, 
alignments of the extended emission-line gas with the radio structures were also found (see iBaum fc Heckmanlll989l ). One 
possible explanation for these observations is that if it is assumed that the EELRs could be bright, large and unresolved from 
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1 they are also known as Extended Narrow-Line Region (ENLR) when they are referred to Seyfert galaxies. 
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the host galaxy, the optical shape of the object would be dominated by the EELR in the case of optically poor spatially 
resolved radio galaxies. So, this is a strong indication that the radio-jet could be the triggering mechanism for the ionisation 
of these regions. 

With the higher resolution instruments, several studies were carried out with the aim of understanding what kind of 
mechanisms can modify the physical condition of the kpc-scale interstellar or intergalactic medium (ISM or IGM) . Specifically, 
any ionizing mechanism shoud be able to explain the emission of radiation in a place with no apparent local source of ionising 
photons, where the stellar radiation could not account for the large amount of photons that are being observed. A very wide 
sample of objects, among them Seyfert and radio galaxies, have been imaged by HST and many of them show large amounts 
of highly ionized and disturbed gas on their outskirts. 

Two important issues concerning these extended ionised regions remain uncertain and they are being discussed: their 
origin and the mechanisms involved in their ioni zation. We have started t o investigate the latte r using the most acceptable 
theories to date: photo i onisation, mixed mat ter ( Binette et al. 19961 . 1997 ) and shock-ionisation ( Dopita fc Sutherland! 1995 , 
19961 : Ullen et alj|2008l : iGroves fc Afienlboiol ). 



Photoionisation assumes that the UV photon flux produced by the AGN ionises the line-emitting gas; it is described 
by the dimensionless ionising parameter U, which accounts for ra diation field dilutio n as t he distance from the nucleus 
increases. More complex photoionisation models were performed bv lBinette et al. I { 19961 . 1997 ) who introduced the presence 
of gas clouds with different optical depths, and different values for U. On the other hand, shock-ionisation is triggered by 
shock wa ves driven by the radio jets. Those shock wa ves compress, heat , and accelerate the IGM, t hus causing it to produce 



emission. 



Dopita fc Sutherland! (l995l . ll996l ) and later Ullen et all |200sl ) and lGroves fc Alien! (|201ol ) developed MHD models 



for shock-ionisation, for both low and high shock velocities. When high speeds are considered, a radiative precursor region 
is created in front of the shock wave by the most energetic photons. As a consequence, the emission region gets bigger (it is 
composed by the precursor and the hot shocked gas) and the spectrum shows a mixture of io ni sation and e xcit ation states. 

Both scenarios have been tested in a wide sample of active galaxies. ICapetti et al. | (|l995al rbl, [l996l . ll997l ) and lWinge eTal] 

1 19971 ) were the first to show that shock-ionisation best explains the optical e mission in the NLR in nearby Seyfert galaxies 
(Mrk 3, Mrk 6, Mrk 573, NGC 1068, NGC 4151 and NGC 7319). Recently. [Rosario et all l|2010al |b|) have shown th a t this 
mech a nism also occurs in NGC 5929 and SPSS J1517+3353. In the case of powerful 3CR radio sources, [fiinstein et al. d 19991 , 
2002 ):lMarauez et all j2000h:ISol6rzano-Inarrea et al] |200 j , l2002l ): ls"ol6rzano-Inarrea fc Tadhunterl |2003l ): lTilak et all (|2005l ): 



Christensen et al 



(200C 



Hardcaste et al.l 



201o| ) sIuav. ,[ |- hat tins interaction mil !>«• (.mud in the EELRs of :',( 1 2!H). '■>('■ 2 L 1.1. 



3C 17 1 3C 34, 3C 330, 3C 352 3C 4 35A, 3C 265, 3C 196, 3C 277.3, and also in PKS 2250-41. In contrast. iRobinson et al 



(2000); ISolorzano-Inarrea et al] (|2004j) have proved that 3C 321 and the powerful radio galaxies 0850-206 and 1303+091 are 
photoionised by the AGN. 



1.1 The galaxy 3C 381 

3C 381 is a FR II double-lobed radio source. It has a flux density of 16.6 Jy at 178 MHz (Laing, Riley & Longairlll983T ). 
The radio structures are located in an approximate North-South direction (see figs. 21-22 from lLeahv fc Perlevlll99ll ). the 
radio-axis being defined by the line that connects the radio core with the two hotspot s. 3C 381 is a large ~ 74 arcsec radio 
source, the south lobe is larger than the north one (40.5 and 33.4 arcsec, respectively iLeahv fc Perlevl ll991): its inner edge 
is more separated from the central source than the north lobe is. None the less, a fairly symmetric rotational pattern can be 
distinguished near the central despite the separation. 



The galaxy was studied at optical wavelengths by iMcCarthv. Spinrad fc van Breugell (|1995l ). who showed the presence 
of extended [O iii]A5007 emission. They characterized the object as one with the largest high-ionisation emission- line regions 
amongst the low-redshift (z=0.1605) 3CR sources, which was later confirmed by HST/WFPC2 data (Fig.[T|). The optical and 
radio features are not aligned. However, near the centre both radio lobes are rotationally symmetric distorted and misaligned 
with respect to the radio axis (defined by the passage of the jet) in the same way as the [O iii]A5007 e xtended emission is . 
Moreover, the north-west optical elongation reaches the south extreme of the north lobe (see fig. 54 from Privon et al.l 2008 L 



but t here is also an additional much more extended [O iii]A5007 filament 25 arcsec to the north (fig. 5 from lMcCarthv et al 



1995) that seems to be either embedded in or superimposed on the north lobe, in a direction approximately coincidental with 



the radio axis. The g alaxy was classified as a High Excitation Galaxy (HEG) based on its strong nuclear spectrum and its [O 



m]A5007 luminosity (jButtiglione et alJbood . l20icl ) 



As all galaxies in the 3CRR catalogue, 3C 381 is located away from the Galactic plane ( Laing et al.lll983l ). The extinction 
is low enough (E(b-v) = 0.053 mag; NASA/IPAC Extragalactic Database) that the Galactic reddening correction is negligible. 



2 See the DRAGNs Web site, http://www.jb.man.ac.uk/atlas/, which is maintained at the Jodrell Bank Observatory by J. P. Leahy, 
A.H. Bridle, and R.G. Strom. 
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Throughout this paper, we assume Ho = 73 km s _1 Mpc" 1 , Qmatter = 0.27, ft va cuum = 0.73. Hence, at a redshift of z = 0.1605, 
3C 381 is at a distance of 768 Mpc with a projected linear scale of 3.2 kpc arcsec - . 

We are going to investigate the emission coming from the EELR of 3C 381 by using long-slit GMOS/Gemini spectroscopy 
with the aim of determining which of the ionising mechanisms mentioned above best explain the ionisation state of this 
large-scale gas. 



2 OBSERVATIONS 

2.1 HST/WFPC2 imaging 

3C 381 radio galaxy was observed with HST/MVFPG2 as part of the 3CR Snapshot Survey (PI: Sparks). The log of observations 
is summarized on Table [1] Two broad-band F702W images were taken in 1995 February (140s each one) which account for 
continuum emission from the host galaxy. Other two narrow-band ramp filter FR533N33 (300s) images were taken in 1995 
August. This filter was employed to find ionised gas surrounding the galaxy. The WFPC2 scales are 0.0455 arcsec px 

(146 pc px -1 ) for PC mode, and 0.0996 arcsec px -1 (319 pc px -1 ) for WF mode. 

The F702W filter is close to Cousins's R; these images were reduced as described in lMartel et al. I jl999h and then they 



were registered and added in order to produce one final image with higher signal-to- noise ratio (S/N). Data were reduced using 
the STSDAS package within iraf, which includes the standard WFPC2 pipeline processing and cosmic ray removal. Given the 
redshift of the galaxy and its position on the CCD, the ramp filter produces an [O in]A5007-dominated image (Fig. [TJ right 
panel). These structures are located at PA 155°, they extend up to 2.2 arcsec to the North- West and South-East from the 
nucleus respectively, and they constitute the targets for our spectroscopic analysis. 

As the HST images are flux calibrated data, photometry was done in order to calibrate the exposure time needed to have 
a good S/N ratio for the spectroscopy. This procedure was performed by measuring the [O in]A5007 emission and then its flux 
was scaled to the oxygen line in a quasar-like spectrum, which is the most similar template available. Both from the images 
and this analysis, it is clear that the nebular structures are very extended and they are bright enough to stand out over the 
stellar luminosity. However, our main interest is focused on the regions where the stellar component already vanishes, or it is 
so faint that there is no measurable contribution in the spectra. 



2.2 GMOS spectroscopy 

The long-slit spectrum was obtained in 2005 April using Gemini's Multi-Object Spectrograph (GMOS) at the 8.1 m Gemini 
North telescope as part of the program GN-2005A-Q-37 (PI: Feinstein; see Table[2]for the log of observations) . The instrument 
was set up with a 0.5 arcsec width slit, and the B600-G5303 grating R~1700 centred at about 4500 A, which yields a resolution 
of 0.9 Apx" 1 . Only one exposure of 2400s was taken with 2x2 binned CCD, implying a spatial resolution of 0.1454 arcsec px 
The slit was placed at a position angle (PA) of 155°, along the optical- UV elongation. 

Data were reduced by using the GEMINI-GMOS package reduction tasks within iraf. We followed the usual steps of bias 
subtraction, flat-field correction, wavelength calibration, and sky subtraction. We performed cosmic ray rejection from the 
2-dimensional spect rum by using two different tasks: gscrrej from Gemini package and the Laplacian Cosmic Ray Identification 
( van Dokkumll200lr ). but we found some loss of information after the tasks were (separately) carried out. For these reasons, 
we decided not to apply this correction to our spectrum. However, the presence of a cosmic ray hit over a given profile in the 
extracted 1-dimensional spectrum is clearly detected because of its shape, intensity, and width relative to the line itself, so 
cosmic ray contamination was easily removed from our measurements. 

The spectrum covers ~2700 A in wavelength, from 3800 A to 5800 A. It includes emission-lines of several species, from 
the bluest [Ne v]A3424 to the reddest [O iii]A5007, the latter being the most important feature in the spectrum (Fig. [2]). 
The set of identified lines is listed in Table [3] together with line fluxes relative to that of H/3 for each selected position. All 
the locations are referred to the galactic centre, which was morphologically determined both for the spectra as for the HST 
imaging. Since 3C 381 is an elliptical galaxy, and it is well fitted by a De Vaucouleurs profile, the centre can be determined 
very accurately for the HST continuum image. Concerning Gemini spectrum, the centre was located as the maximum intensity 
point within the stellar continuum. This point matches the maximum of the De Vaucouleurs profile. 

Since Gemini spectrum has a resolution of 0.1454 arcsec px" 1 , it gives a linear scale of about 4654 pc px" 1 , without 
projection correction. Fig. [2] covers the range from H,9 to [O iii]A5007; the intensities of these lines are shown in gray scale; 
the contours have been plotted with the aim of emphasizing the line's intrinsic shapes. 
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3 RESULTS 

3.1 HST/WFPC2 imaging 

The galaxy shows a completely different morphology at broad band compared with the narrow band. While the first shows the 
stellar component dominating the flux, the narrow band filters are sensitive to the ionised gas. Its elliptical nature is clearly 
seen in broad-band filter (Fig. [TJ left panel). The elliptical structure is homogeneous and it is well fitted by a De Vaucouleurs 
(r 1//4 ) surface brightness profile. It vanishes below the noise in the HST images at 1.7 arcsec from the nucleus, the stellar 
component being negligible as from this point. 

On the other hand, the HST/WFPC2 ramp filter image shows the nebular emission dominated by [O m]A5007. These 
structures are elongated and they vanish in the outskirts of the host galaxy, giving a total extension from one end to the 
other of about 4.5 arcsec at a position angle of 155° (counterclockwise from the North). The shapes of both features are quite 
different. The NW EELR is narrower than that of the SE; it has its maximal emission concentrated over its inner region from 
where it adopts an S shape towards its outer limit, reaching the south extreme of the northern radio lobe. The SE EELR, in 
contrast, is clearly detached from the central object. The peak of its emission is located at 1.1 arcsec from the nucleus (Fig.[TJ 
right panel). 



3.2 GMOS spectroscopy 

The high-excitation nature of 3C 381 established on the basis of its nuclear spectra is emphasized with our data, where high- 
and low-excitation lines coexist in a region that extends up to 3.2 arcsec (10 kpc) from the nucleus in each direction. However, 
the [O m]A5007 emission is very intense, so it is measurable at a larger distance from the nucleus with a higher S/N than any 
other species. It extends up to 4.8 arcsec to the NW and 7.3 arcsec to the SE. 

The angular size of the entire EELR measured by the [O ill] A5007 line reaches 12 arcsec (38.4 kpc), an extension of almost 
three times what the HST had detected, which indicates the existence of very low surface brightness structures towards the 
outer edges of the EELR. 

The line shapes in the long-slit spectrum show a strong core over the central region, with a blueshifted wing to the NW 
and a redshifted wing to the SE (Fig. Given the complexity of the lines, we have extracted one 1-dimensional spectrum 
for each spatial pixel across the slit to avoid the loss of spatial information. Once the extraction was performed, we found a 
complex spectral structure, where the line profiles can be decomposed into two or even three (in the case of [O in]AA4959, 5007) 
Gaussian components. 

The Gaussian decomposition was performed by using the ngaussfit task from the STSDAS package and the splot task from 
the NOAO package within iraf. The total flux together with its rms error were obtained with ngaussfit and then checked 
with splot. The ngaussfit task works by iteration based upon a set of initial conditions on its three free parameters: central 
wavelength, intensity at the peak and fwhm for each component, in addition to the baseline continuum. We set initial 
conditions for the central wavelengths and line intensities; the fwhms were left to converge by themselves. Finally, the initial 
values we proposed did not change significantly over the iterations and all the profiles reached their best fit with narrow 
(fwhm < 800 km s _1 ) components for both forbidden and permitted lines. The line intensities decrease from the nucleus to 
the outskirts for the north-west EELR, but there is a sudden increase towards the south-eastern region at around 1 arcsec 
(~ 3.2 kpc). This behaviour was detected in the most intense lines as well as in the weakest ones. It reveals that the emission 
increase towards the SE in the [O m]A5007 contour-map (Fig. [1] right panel) is common to the whole EELR SE-gas. 

Strong line-splitting was found in oxygen forbidden lines but also in hydrogen recombination lines. Although 3C 381 has 
been considered as a Broad-Line Radio Galaxy (BLRG) on the basis of a weak broad component in Ha ijGrandi &: Osterbrockl 
1978), no broad component was needed to reconstruct the profile in any Hydrogen line. Even when Ha lies outside the 



wavelength range we cover, we are able to sa y that 3C 381 has a Seyfert H-like sp ectrum. It has been also confirmed not only 



by X-ray nuclear data taken with Chandra (jHardcastle. Evans fc Croston 200S, and reference s therein), but also by recent 



optical nuclear data taken with Telescopio Nazionale Galileo (TNG) ( Buttiglione et al. 20091 ') . which shows no evidence of 
such a broad component. For these reasons we assume 3C 381 is a Narrow Line Radio Galaxy (NLRG), which is in keeping 
with its lack of a quasar-like core. 

In addition to this line-splitting, the strongly distorted [O m]A5007 emission- line profile also has a broad underlying 
component, mainly manifested inside the inner ±3 arcsec (±10 kpc). The broad wings in [O in]A5007 (which are also detected 
in [O in]A4959) make it necessary to use more than a single narrow gaussian fit. The narrow component can be easily followed 
throughout the spectra, and its parameters are measured with high accuracy. On the other hand, the errors associated with 
the broad component parameters are high, specially the fwhm. The measurement of this component is much more difficult 
than that of the narrower one, but its underlying presence is not arguable. 
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3.2.1 Kinematics of the emitting gas 



Due to its brightness, the [O m]A5007 emission- line was used to compute the velocity field relative to the host galaxy systemic 
velocity (Fig. [3]). A fairly symmetric distribution to each side of the nucleus is distinguished as an organized and collective 
pattern, coming from a gas component probably related to a disk. This rotation curve tends to stabilize around 300 km s _1 
at ~1.5 arcsec (~ 5 kpc) from the g alactic centre which is consistent with gravitational motion within the galactic potential 



i Tadhunter. Fosburv fe Quinnl 1 19891 1. However, there are also high velocity components superimposed on that curve which 
follow neither ordered nor collective behaviour. Such a velocity shift cannot be explained under the gravitational hypothesis. 
The highest velocity components are located within the inner ±3 arcsec (±10 kpc); they are separated from the rotation curve 
by ±300 km s _1 approximately. They deserve special attention because those disturbed kinematic components are also the 
broader ones (Fig. 0). 

The high-velocity displacements among the kinematic components could be related to disturbing processes triggered by 
the radio-emitting components in different scenarios. Evidence of such processes has been found not only in radio galaxies 
where the radio components o verlap the line-emiting gas, bu t also in radio sources where the radio structures are found far 
beyond the innermost EELR (|Sol6rzano-Inarrea et al.l 120011 ). The same happens in 3C 381; the radio and optical features 
have different scales, excepting the outer edge of the NW EELR that reaches the s outh extreme of the no rthern lobe, in 
addition to the other ionised filament, detected in the same location of the north lobe iMcCarthv et al.l (|1995). 3C 381 is one 
of these galaxies where the interaction between radio and optical features can be identified even when these components are 
not overlapped. 

The kinematic disturbed components we find are strong indicators that this gas is undergoing a violent disturbing process, 
very probably triggered by the radio structures. So, any modelling of the undergoing emission process must also account for 
the energy that is accelerating this component of gas. 



4 DISCUSSION 



To understand and identify the main ionising mechanism of the EELR, we test the line-ratios from the 3C 381 spectrum 
with the results predicted from modelling of several mechanisms that could ionise the interstellar gas. In this process we have 
to identify the main source of energy of the line-emission, which could be the AGN (e.g. AGN-photoionisation, mixed-matter 
models) or a local phenomena (eg. shock- ionisation). Most of these mechanisms are already modeled with predicted line-ratios 
that can be easily contrasted with the observed data. 

A direct way to model the proble m of the ionising ph oton budget is to calculate the dimensionless ionisation parameter 



U through the empirical relationship ( Penston et al. 1990h over the entire nebular region 



U = -2.74 - log 



[Oil] A3727 \ 
[OIII] A5007 ) 



This is the simplest way to understand if the UV photons coming from the AGN could be responsible for these ionised 
regions of extended gas. This relation was obtained under the assumption of solar abundances and gas density of 100 cm -3 . It 
covers all possible configurations applicable to AGN, with spectral indices a (<!>„ ~ v +a ) varying in the range of —2 < a < — 1 
and black body temperatures in the range of 10 s K < T < 2.1 x 10 7 K. The ionisation parameter is used to find an estimation 
for the amount of ionising photons emitted per time un it by the AGN (Qo), for which the ambient gas has to be known 
(Qo = U4ivr 2 cn,,). iRobinson. Tadhunter fc Dyson |2002l) performed long-slit spectroscopy and determined that at 1.4 kpc 
from the nucleus (PA=155°) the density has a value of 370 cm~ 3 . 

We measured Qo by using this local density together with U calculated in the same location where the density was 
obtained. At this point (1.4 kpc), U = 9.35 x 10 -2 . As our long-slit spectrum allows the spatial measurement of [O n]A3727 
and [O iii]A5007, we can compute the line-ratio and, under the assumptions of Penston's relationship, trace the spatial 
evolution of U (Fig. 4). Even if the density changes as a function of position, Qo is supposed to be constant. So, we obtain 
Qo = 2.4 x 10 56 photons s _1 . Both Qo and U must be considered as a minimum because the inclination is unknown, and it 
has not been corrected. 

This Qo is c omparable to that of 3 C 273, the most powerf ul radio-loud quasar known at low redshift (Qo = 3.8 x 



smparable to that or 30 Hi, the most p owertul radio-loud quasar known at low redshht (Qo = 3.8 x 
Robi nson et al. l (|200d i% lHardcastle et all (|2009l ) have shown that what they called X-ray 'accretion related' 



10 SB photons s" 

nuclear luminosity (absorption corrected) is a good indicator of the AGN power. The X -ray luminosity of 3C 381 is almost an 
order of magnitude lower that that of quasars (see fig. 13 from Hardcastle et al. 2009) (and the [O iii]A5007 luminosity too). 
So it is very likely that the Qo previously calculated is a very overestimated measurement of the ionizing power in 3C 381. 



In addition, other problem arises concerning U. We expect that U evolves as U ~ n e 



So it would be sensible to 



3 their published value is Qo = 3 x 10 55 photons s 1 ster 1 
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the geometrical dilution of the central ionising radiation field, and the changes in the density. Fi g. [4] reflects the behaviour of 



U, through [O m]A5007/[O n]A3727, under the assumptions of constant density ( Penston et al. 1990h . It can be clearly seen 
that U does not follow the r~ 2 dilution over a large distance. In contrast, U increases towards the south-eastern EELR and, 
even more, we must highlight that every emission- line in the spectra increases its intensity towards the SE (at ~ 1 arcsec). In 
order to produce such an increase in the emission (as U shows), the density should decrease faster than r~ 2 over a kpc-scale 
distance, which seems to be unlikely. It is worth noting that the [O in]A5007/[O n]A3727 line-ratio will significantly change 
its predicted behaviour (with respect to the single ionising source) if there is other off- nuclear source of ionising photons, such 
as the local radiation field created by local radiative shocks. We will discuss it later. 

Therefore, from the overestimated ionising budget and the observed behaviour of the [O m]A5007 and [O n]A3727 line- 
ratio, we think that it is very unlikely that a simple AGN-photoionisation could be the dominan t ionising mechanism . 



Then, a more complex photoionisation model is required. A clever approach was developed bv lBinette et al. (1996) doing a 



more s ophisticated modelling of t he way in which a clumpy medium absorbs the impining radiation field. In this mixed-matter 
model ( Binette et al. 19961 . 1997 ). the EELR is composed by two different cloudy systems. The matter-bounded (MB) clouds, 



which are the optically-thin lower-density cloudlets, absorb a fraction of the incoming radiation. Since the highest-energy 
photons are absorbed here, this is the place where the highest-excitation lines are formed. The lonisation bounded (IB) clouds 
absorb the diluted ionising field that exits the first component. The lowest-excitation lines in the spectrum are formed in this 
optically-thick cloudlets. In this clumpy system, the radiation field emitted by the central surce is partially absorbed by the 
MB clouds, and once it has been leaked, it is completely absorbed by the I B clouds due to their higher density. The extension 
of the MB component is defined by the avail ability of He + ioni sing photons ( Binette et al. 19961 ). which in turn are manifested 
by the presence of Hen recombination lines (|Osterbrock 1989h . 



However, the He nA4686 emission 3C 381 extends up to 2.5 arcsec (~ 8 kpc) over each EELR, without the strong variation 
that the model proposes. The line intensity at the end of each EELR drops to ~ 10% of its maximum value at the galactic 
centre, while the intensity of [O m]A5007, the main feature in the spectrum, drops to less than ~ 5% of its maximum value in 
the same spatial extent. In fact, the extent to which the line is measured is in itself a problem, since it sets the spatial limit for 
the MB clouds. Consequently, there is no mixed-material inside the inner 2.5 arcsec (~ 8 kpc); therefore, there applies again 
the hypothesis proposed by the simple AGN-photoionization model discussed before, which fails in reproducing the scenario 
in 3C 381's environment. 

The modelling is implemented by diagnostic diagrams parameterised by Am/i (defined as the solid angle subtended by 
the MB clouds with respect to that of IB clouds) which varies as the He n/H/3 line-ratio does. The diagram plotted in Fig. [5] 
uses the oxygen ratio [O ll]A3727/[0 in]A5007 against Hell/ H/3, which became an excitation axis given its relation with A M /j. 
It is worthwhile comparing this diagram with fig. 7 from Binette ct al. I |l996h to ensure that the observed trend actually 



responds to changes in U rather than Am/i (see also iFerruit et aL 19991 . for the same kind of analysis on Mrk 573). Each 



EELR was plotted separately from the other, connecting the measurements to follow their trend with respect to the A M /i 
sequence. The NW measurements were plotted as dots, and those from the SE as triangles; the central spectrum is shown by 
a square. 

It is very unlikely that photoionisation models could account for the whole observations. Even when there should be a 
contribution from photoionisation, neither the observed increment in U toward the SE nor the line-ratios can be reproduced 
by these kind of models. 

A plausible explanation for an increasing ionisation parameter (note that U increases because [O ll]A3727 increases as 
well as [O in]A5007 does) is the shock- ionisation model, since either a local source of ionising photons is needed or the 
AGN has improved its photoionising power as a consequence of the disturbing processes undergone by the EELR material 
I Solorzano-Inarrea fc Tadhunter 2003). This mechanism, also known as jet-cloud interaction, might provide a local ionising 



field if the expanding shock waves are fast enough, but such a scenario should be manifested not only in line-intensities but 
also in the kinematics of the region. 

When a shock wave is capable of ionising the medium where it propagates, the emission of radiation usually comes from 
the recombination region behind the shock front. But the fastest shocks (v > 170 km s _1 ) create UV photons that can diffuse 
the jet both upstream and downstream, changing the ionisation conditions of the entire region. Photons travelling upstream 
create an Hn-like region (known as the precursor) ahead the shock, but also the post-shock gas is affected by photons travelling 
downstream. As a consequence, the line-emitting region gets an i mportant contrib ution from the precursor: it is 10% to 20% 
brighter than the shock in H/3 (jPopita fc Sutherland!! 19951 . Il9961 : lAllen et allboosh . 

These models are also implemented through line-ratio diagnostic diagrams. We have drawn (Figures [6] to 1 1 1 j) the 
shock +precursor g r ids (d otted lines, one line per magnetic paramete]0) from the most recent modelling of these phenom- 
ena by Allen et al. 1 2008T ). Models with solar abundance and pre-shock gas density of 100 cm -3 were chosen. The only-shock 



4 The magnetic paramete r is defined as B/n 1 / 2 , B is the transverse magnetic field, and n is the pre-shock gas density 
jDopita fc Sutherland1ll995l) . 
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grids (no precursor) lie ou tside the range of our observations (compare the location of our data in Fig. [6] with fig. 1 from 



Dopita fc Sutherland 19951 . where both grids were plotted together). Shock velocities are also shown (dashed lines). Data are 



plotted as explained before. 

We have used both [O m]A5007/H/3 and He n/H/3 as excitation axis. These ratios, as well as the other pairs of lines 
employed, are close in wavelength to minimise the unknown reddening effects. In other c ases we have used lines from the sam e 
species, to reduce the consequences that an unknown metallicity may have in our results jBaldwin. Phillips fc Terlevicb|[l98ll 



As a general characteristic, our [O lll]A5007/H/3 and He ll/H/3 values are higher than predicted by the model. In addition, 
our [O m]A5007/H/3 measurements are also much more constant. Anyway, the fit is good despite these minor inconsistencies, 
and the best fit is reached in Fig. [7] which employs those line-ratios. 

The highly excited [Ne v]A3424 emission- line is used toghether with [O ll]A3727 in the next two plots (Fig. [8] and [9}. 
The [Ne v]A3424 emission extends up to ~ 5.5 kpc away from the AGN; since ~ 95 eV photons a re needed to produce th is 
emission-line, it is very difficult for any photoionisation model to explain such an extended emission jKomossa fc Schuljl 19971 ). 
In contrast, jet-cloud interaction may provide local energetic photons. As it is seen from the plots, these models also success 
in reproducing the relative intensity of [Ne v]A3424 with respect to the lower excited [O ll]A3727 emission- line. 

The last two diagrams (Fig. [TOl and lll|) use the temperature-sensitive line-ratio [O in]A4363/[0 m]A5007, which are 
particularly complicated to be reproduced by any kind of model. Since the "temperature problem" associated with this ratio 
( Dopita fc Sutherla nd 1995) is widely known, we do not draw any conclusion based upon these plots, but we show them 
because, despite the fact of having poorer fits, the data behave as in previous diagrams. 

Both high- and low-excitation emission lines are well fitted by the shock-ionisation model, although there is a clear 
underprediction regarding to [O in]A5007 and He II. Furthermore, the extent to which the highest excited line in the spectrum 
is measured (i.e. [Ne VA3424]) requires a budget of energetic photons that can be locally created by this kind of fast radiative 
shocks. In fact, the behaviour of the ionisation parameter U is also explained by a local ionising field. Data tend to be located 
over the fastest (500 km s _1 < v <1000 km s _1 ) shocks with no exception. It is worth highlighting that the disturbed [O 
iii]A5007 components reported on the velocity field move with velocities in the same range as those predicted by the shock- 
ionisation model. This is a strongly remarkable point since it is independent evidence to support the shock-ionisation idea as 
the main source of the ionisation of the gas. 

Certainly the most spectacular interaction between the radio jets and the IGM is observed at the edges of the lobes, 
where the shock waves driven by the jet give raise to the bright, highly overpressured hotspots. The effects of shock waves 
in the way toward the hotspots have been also det ected in the optical-em itting gas. An impressive example of both jet-cloud 
and lobe-cloud interaction was reported in 3C 299 jFeinstein et alll999l ). where the EELR overlaps the radio lobe. However, 
a turning point was then achieved bv lSolorzano-Inarrea et al* I 200ll ) who showed that, despite the large distance in-between 
the EELR and the radio lobes in their sample, the EELR's gas undergoes jet-induced shock-related disturbing processes. In 
these latter cases the key evidence came from the strongly disturbed kinematical patterns observed within the gas. The same 
kind of disturbing motion that we find in 3C 381. 

Jet-driven shock waves are triggered near the AGN (near compared with the lobe-scale) since the jet is supersonic in 
these regions (Leahy; www.jb.man.ac.uk/atlas/dragns.html). Conse quently, those shocks g ive raise to high-pressured regions 



within the clumpy EELR. Such a pressure-excess was reported by iRobinson et al.l (|2002l ). They made use of their density 



measurements (we have used it at the beginning of this Section) to show that the EELR's gas is overpressured with respect 
to the X-ray halo. They claimed that jet-induced shocks might be playing an important role in the region, but no kinematical 
argument was then available to support this scenario. The velocity field we have shown in Section fa. 2. II accounts for such an 
argument. 

We have shown that the jet-driven shock-ionising waves that are able to reproduce the observed line-ratio are also those 
that move with velocities entirely compatible with the high-velocity line-splitting shown by the velocity field. So, the kinematics 
of the region, its ionisation and also its pressure excess are now explained as consequences of the same process: the interaction 
of the radio-jet with the IGM. 



5 CONCLUSIONS 

From the study of the long-slit Gemini/GMOS spectra over the extended emission-line region found in the HST imaging of 
3C 381, and the state-of-the-art modelling for the energetic balance for these structures, we find the following results: 

- The nebular region extends up to 12 arcsec (~ 38 kpc) and it is larger than the stellar component as seen in the HST 
data, which vanishes below the noise at 1.7 arcsec (~ 6 kpc) from the nucleus. 

- The kinematics of the EELR shows features of disturbing processes, where high velocity gas components are clearly 
found not only near the central region but also towards the outskirts of the EELR. Differences up to 600 km s _1 in velocity 
are observed within the nebular structure. 

- Line-ratios have been used to compare theoretical predictions with our data and we found that the shock-ionisation 
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model with a precursor component predicts the observed line ratios. The fastest shocks (500 km s _1 < v < 1000 km s _1 ) 
are capable of reaching the extreme values of our data; there are variations amongst diagrams, but there is no contradiction 
between any line ratio. 

- Shock-ionisation models are likely to be the dominant ionising process taking place within the EELR. In addition to the 
diagnostic diagrams that reproduce the observed line-intensities with shock velocities coincident with what the velocity field 
shows, the local ionising photons created by the shocks themselves provide a plausible explanation for both the existence of 
very extended high- excitation lines and the increment of U despite the increasing distance. 

This paper has been typeset from a TeX/ DTeX file prepared by the author. 
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Table 1. Log of HST observations 



Camera 


Date Exp. Time 


Filter Name 


Exposures 






WFPC2 1995 Feb 11 140s 


F702W 


2 






1995 Aug 13 300s 


FR533N33 


2 






Table 2. Log of Gemini observation 










Instrument 


Slit Grating 


Date 


Exp. Time 






CMOS 


0.5 arcsec B600.G5303 


2005 Apr 7 


2400s 






Table 3. Set of identified emission-lines : 


in 3C 381 spectrum. 






Line^ 


Flux 3.2 arcsec NW 


Flux 1.6 arcsec 


NW Flux GC^ 


Flux 1.6 arcsec SE 


Flux 3.2 arcsec SE 


[Nc ni]A3345 




0.18 


0.47 


- 


- 


[Ne v]A3424 




0.81 


1.05 


0.57 


- 


[O n]A3727 


9.93 


4.28 


2.58 


3.01 


4.91 


[Nc ni]A3869 




1.47 


1.19 


1.16 




H<5 (4100) 




0.22 


0.19 


0.19 




H 7 (4340) 




0.54 


0.47 


0.52 




[O m]A4363 




0.27 


0.38 


0.29 




He n (4686) 




0.26 


0.31 


0.29 




[Ar iv]A471l£, 




0.04 




0.11 




[Ar iv]A4740~" 






0.07 






H/3 (4861) 


1 


1 


1 


1 


1 


[O m]A4959 


4.24 


4.1 


4.36 


4.53 


4.33 


[O m]A5007 


13.86 


12.42 


13.32 


13.57 


12.57 



a Wavelengths in Angstroms 
" Galactic centre 

c Possible contamination from [Ne iv]A4718 
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Figure 1. HST/WFPC2 images of 3C 381. Left: The stellar component accounting for the elliptical structure. Broad band filter image. 
Right: [O m]A5007 -dominated image shown in grey-scale and contour map. The contours levels a re 5%, 7%, 10%, 18% and 22% of the 
maximun intensity at the nucleus position. The radio axis direction (PA=4°, dc Koff ct al. 1996) is also shown at the bottom right of 
this panel. 
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Figure 2. Two-dimensional Gemini spectrum, from H/3 to [Om]A5007. The ordinate scale indicates the SE (negative) and NW (positive). 
The contour levels are shown only to enphasise the grey-scale contrast. 
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Figure 3. Velocity field computed from the components that form the [O iii]A5007 emission-line profile. As in Fig. [2] SE is plotted as 
negative coordinates and NW as positive. 
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9 f Z 

^2^gy[iio]/^ooqy[iiio] 

Figure 4. Spatial variation of the ionisation parameter U expressed through the oxygen ratio iPenston et al.iri99Clh .As in Fig. [2] SE is - 
and NW is +. Note that U does not follow the r~ 2 geometrical dilution but grows up at ~1 arcsec towards the SE. Every emission-line 
in 3C 381 's spectrum increases its intensity at the same location. 
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Figure 5. Continuous line: mixed medium photoionisation model from lBinette et"al The black cross shows the position where 



i-M/I 



1. Its extreme values are shown at the end of the sequence. Each EELR was plotted separately, the data are connected to follow 



their trend with respect to the A M /j sequence. The NW measurements were plotted as dots, and those from the SE as triangles; the 
central spectrum is shown by a square. 
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Figure 6. Dotted lines: shock-ionisation + precursor models from I Allen et all ll2008f ), one line per magnetic parameter (B/n 1 / 2 ) value. 
Shock velocities are drawn as dashed lines. See Fig. [5] for data reference. 
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Figure 7. This plot shows the relation between two excitation axes, [O in] A5007/H/3, and the other one implemented bv lBinette et all 
jl996h He ll/H/3. See Fig. [6] for model references. 



© RAS, MNRAS 000, [T]-?? 



18 V. Reynaldi and C. Feinstein 



~\ 1 1 1 1 1 1 1 1 1 r 



300 km s 



200 km s" 



CD 
CO 

> 
0) 



cv 
co 



CD 
o 



o 



2 



_l , , 400 km s" 



500 km s" 



700 km s 



.<••/ l\ \ 

7 ', / 
; f 
\ / 
V 



1000 km s 



-l 



B/n 1/2 = 0.5, 1, 2, 4, 5, 10/xG cm 



3/2 



_l I I I L_ 



1.5 



0.5 







0.5 



log(HeIIA4686/H /? ) 



Figure 8. The high-excitation lines [Nc v]A3424 and He II are combined with the lower excited [O n]A3727 and H/3. The shock-ionisation 
models can reproduce the observations when the highest velocities are considered. 
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Figure 9. The same than Fig.[8]with [O ni]A5007/H/3 as the excitation axis. See Fig.[6]for model references. 
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Figure 10. The temperature-sensitive line ratio [Om] A4363/[Oin] A5007 against [Oin] A5007/H/3. Despite the temperature problem, the 
diagram shares the same characteristic of the previous plots. See Fig. [6] for model references. 
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Figure 11. The temperature-sensitive line ratio [O m]A4363/[0 m]A5007 against He ll/H/3 as excitation axis. See Fig 
references. 
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